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Weakly disordered two-dimensional superconductors undergo a Kosterlitz-Thouless (KT) transition, 
where at a critical temperature vortices proliferate through the system and destroy the supercon- 
ducting (SC) order. On the other hand, it was suggested that for large disorder the systems separates 
into regions of high SC order, and it is the percolation of coherence between these regions that is lost 
at the critical temperature. Here we demonstrate that these two descriptions are just the dual of 
each other. A vortex causes loss of local correlations, and thus the loss of percolation of correlations 
is concomitant with percolation of vortices on the dual lattice, in the perpendicular direction, i.e. 
the KT transition. 

PACS numbers: 74.20.-z, 74.45. +c,74.81.-g 



The interplay of disorder and superconductivity has been 
a subject of research from the early days of BCS the- 
ory |l| , when Anderson has demonstrated that weak dis- 
order does not affect the BCS critical temperature 0- 
It was later argued that superconductivity indeed per- 
sists even when the single-particle states are localized 
but eventually, with strong enough disorder, super- 
conductivity is destroyed. Thus the critical tempera- 
ture is reduced with increasing disorder, until it is sup- 
pressed all the way to zero for large enough disorder, 
indicating a zero-temperature transition from a super- 
conducting (SC) to an insulating phase. The situation 
is even more intriguing in two dimensions, where, even 
without disorder, there could be no long-range SC or- 
der at any finite temperature T due to the Mermin- 
Wagner theorem 0. Nevertheless, one expects a finite- 
temperature Kosterlitz-Thouless (KT) transition [5] be- 
tween a phase with power-law decaying correlations to 
a phase with exponentially decaying correlations. This 
transition is driven by the unbinding and proliferation 
of vortex-anti vortex pairs, which destroy the SC order. 
The effect of disorder on the KT transition has become 
even more relevant since the experimental observation 
of superconductor-insulator transition in thin disordered 
films, about two decades ago [1,0] ■ Weak disorder should 
not affect the transition, according to the Harris criterion 
Q. In the presence of strong disorder, however, it has 
been established theoretically and observed exper- 

imentally [Hill that the SC order parameter fluctuates 
strongly across the sample, creating "SC islands", where 
the SC order is high, surrounded by areas of weaker SC 
correlations. Thus, with increasing temperature the co- 
herence between neighboring SC islands is quenched, un- 
til percolation of coherence from one side of the sample 
to the other is lost, leading to the loss of global SC or- 
der dum. This description predicts that local SC order 
may persists even when global SC order is lost, consistent 
with recent experiments [l7l - [l9j ]. 

In this letter we study the thermal phase transi- 



tion from a SC to a normal state in disordered two- 
dimensional superconductors. In particular we address 
the question of relevance of the KT description to this 
transition, and confront this description with the percola- 
tive description. To that end we begin with the negative- 
U Hubbard model, 
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where indicates a sum over nearest neighbors, Cj a 
creates a spin-cr electron at site i, Uj is the hopping in- 
tegral, taken to be the unit of energy in the following, 
and U < is the on-site attractive potential. The site- 
specific disorder Vi is taken from a uniform distribution 
of width 2W such that V l € [-W, W], but for each re- 
alization we ensure that J^- Vi = 0, and the chemical 
potential fi determines the average density n. The choice 
of the Hubbard model is because it can lead, depending 
on parameters, to a BCS transition, to a KT transition or 
to a percolation transition, and thus it is general enough 
not to limit a priori the possible transitions, unlike, e.g. 
the XY model, which is known to give rise to a KT tran- 
sition [B| and may bias the results towards this particular 
transition. 

In order to consider SC fluctuations, crucial for the 
description of the transition, we employ a method that 
takes into account thermal phase fluctuations, but ig- 
nores the quantum ones [la, l2(ij ]. In short, apply- 
ing a Hubbard-Stratonovic transformation to the Hub- 
bard Hamiltonian {!]), with a local complex Hubbard- 
Stratonovic field, Aj, and ignoring the temporal depen- 
dence of these fields (quantum fluctuations), the partition 
function becomes: 
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Z = Tr[e-P H ] = J X>({A i5 A*})Tr/[e 
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with the Bogoliubov-de Gennes Hamiltonian [2l| 
HsdG({Aj}) given by 

H-BdG = ~J2 taClCfr + J^iVi-riClCia 
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Here Try traces the fermionic degrees of freedom over 
the single-body Hamiltonian H BdG and can be evaluated 
exactly using its eigenvalues 22]. The integral over the 
fields {Aj, A*} can then be calculated using the (classi- 
cal) Metropolis Monte-Carlo (MC) technique [23| . One 
should note that unlike the usual BdG approach, here 
A, are auxiliary fields, and except at zero temperature 
where the saddle point evaluation of the partition func- 
tion gives rise to the BdG solution, they are generally 
different from the local SC order parameter < C^C^ >. 

This procedure allows calculation of any ensemble- 
averaged quantity. We first calculate the average cor- 
relation function Dy =< cos(f5j — 6j) >, where 0, is the 
phase of the SC order parameter < C^CVf >, and i and 
j are on the opposite edges of the sample. Bulk SC order 
exist only if these phases are correlated, hence the decay 
of to zero signals the loss of bulk SC order in the 
system. This calculation was carried out in two ways. 

can be calculated directly using the MC procedure. 
Alternatively, we first calculate the whole distribution of 
9i — 9j , for any pair i and j . We then compare this distri- 
bution, for every temperature, to the one expected from 
a single Josephson junction, 
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-Jij(T) cos(0i-0j). 
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The excellent fit to the data (not shown) allows us to 
determine the (temperature dependent) effective Joseph- 
son coupling Jij(T). Once this parameter is determined, 
D^ can be directly evaluated, using the Hamiltonian @ 
to be Ii(Jij/T)/I (Jij/T), where I n (x) is Bessel function 
of the first kind. Both approaches give almost exactly 
the same results, plotted in Fig. [TJ As can be seen in the 
figure, increasing disorder suppresses the critical temper- 
ature T c , as expected. Except for the change in T c , the 
decay of correlations at finite disorder looks very similar 
to the zero disorder case, expected to be described by the 
KT transition. 

In order to quantify this observation further, we cal- 
culate the sample-averaged vortex (and anti- vortex) den- 
sity, n v . The local vortex density for a single snapshot 
during the MC procedure is determined by the clock- 
wise integration of the phase of the SC order parameter 
around a single plaquette. The value of this integral can 
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FIG. 1: The edge-to-edge phase correlation Dij (left axis) and 
the vortex density n v (right axis) as a function of temperature, 
for different values of disorder for 18 x 18 system, and U = —4, 
and values of disorder depicted in the legend. With increasing 
disorder T c , where bulk coherence is lost and vortices start to 
proliferate, shifts to lower values, but the curves still look 
similar. The circles denote the percolation temperature T p 
(see text). 



assume only integer multiples of 27r, indicating the ex- 
istence of vortices (n > 0), anti- vortices (n < 0), or no 
vorticity (n = 0). In Fig.Q]we also plot n v as a function of 
temperature. As can clearly seen, the loss of global phase 
coherence occurs when vortices start to proliferate, for all 
values of disorder, indicating a KT-type transition. 

In order to further identify the KT characteristics of 
the transition, we adopt a particularly elegant analysis 
proposed by Polkovnikov, Altman and Demler [H)], and 
used in experimental analysis by Dalibard et al. [2o| . 
The idea is to perform the double sum over all pairwise 
correlations D^, up to the system size L, 
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where the second relation defines the exponent a. As 
T^Owe expect Dy — > 1 meaning that A(L, T — > 0) = 1 
and therefore a(T — > 0) = 0. In the opposite limit, when 
T > T c then D^ decays exponentially with the distance 
\i - j\, leading to A(L,T > T c ) -> L~ x and a -> 0.5. 
At the KT transition (see (24|), a universal jump from 
a = 0.25 to a = 0.5 is expected. 

Fig. [2] depicts our results for the PAD exponent a as 
function of temperature, for different values of disorder. 
All the curves (except for the highest disorder, proba- 
bly already on the insulating side) exhibit a very similar 
behavior - a moderate rise from a = to a ~ 0.1 and 
then a rather abrupt rise towards a ~ 0.5, demonstrating 
a change from power-law decay of correlations to expo- 
nentially decaying correlations, again consistent with the 
KT description. We attribute the deviation of the value 
of the exponent where the abrupt change occur from the 
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FIG. 2: The PAD exponent 01 a for u = ~ 4 > 18 x 18 system 
for different values of disorder, exhibiting an abrupt change, 
signalling the transition. Large dots correspond to T p from 
the percolation picture (see text). 



expected value (a — 0.25) to the finite size of our sys- 
tem, as the coherence length £ j which is approximately 
given by ftvp/A, where vf is the Fermi energy and A the 
SC gap, is not much smaller than the system size (about 
a factor of 2 or 3). The fact that the same behavior is 
observed at zero disorder, where the KT description is 
indeed expected to hold, supports this explanation. 

Having established the relevance of the KT descrip- 
tion for finite disorder, we now turn to check whether the 
same data can be described by the percolation picture. 
According to the discussion above, we can determine the 
effective Josephson coupling Jy between nearest neigh- 
bors. With increasing temperature Jy decays, until cor- 
relations are lost, around ■Zy(T) ~ T (corresponding to 
Dij ~ 0.45). Thus we can set this value as a threshold 
for an existing coherence link between nearest-neighbors, 
and determine the temperature where edge-to-edge per- 
colation is lost, T p . Fig. [3] demonstrates the connectivity 
of the network for temperatures from below to above the 
percolation temperature, where the existing links are col- 
ored by strength (of Jij or Dij), and the missing links are 
those with Dy below the above threshold. The percola- 
tion temperature T p can be determined for each disorder 
value, and, in fact, corresponds to an average link den- 
sity of 1/2, consistent with bond percolation on a square 
lattice. T p is also depicted as circles in FigJTJ demonstrat- 
ing that it is indeed in the region where bulk correlations 
decay and vortices start to proliferate. 

In order to quantify this last statement further, we 
have also plotted in Fig. [3] the local vortex density n v , 
in brighter shades of green for higher density. As can 
be clearly seen, as the local n v increases, the neighbor- 
ing decrease and disconnect. In Fig. [4] we plot the 
dependence of the local Jij on the local n v , for different 
values of disorder. This dependence is observed to be 
independent of disorder (except for small corrections at 
zero disorder and a slightly wider distribution with in- 
creasing disorder), leading to the local Jij, up to some 




(a) T = 0.02, deep in the SC (b) T = 0.048, in the SC state 
state close to the transition 




(c) T = 0.052 insulating state, (d) T = 0.074 deep in the 
just crossed the transition insulating state 



FIG. 3: Percolation transition for a 18 x 18 system (U = —4) 
at disorder W — 1.3: Nearest neighbor links Dij (lines) col- 
ored from orange for strong correlations, to black for weaker 
correlations (close to cutoff dij = 0.45), and disconnected for 
correlations below the threshold. Temperature changes from 
deep in the SC regime (a) to deep in the insulating regime 
(d). Local vortex density n v (background) colored from dark 
green for low density, to bright yellow for high density, demon- 
strating the correlation between the local Dij and the local 
n v (see also Fig. [4]) . 

small fluctuations, being a unique function of the local 
vortex density. This observation allows us to relate the 
two pictures: in order for the percolation path of the Jy's 
to disconnect, a perpendicular line connecting points of 
n v , larger than some threshold value, has to be formed. 
Thus loss of percolation of the Jij, corresponds, in fact, 
to percolation of the vortices in the dual lattice in the 
perpendicular direction. Since this necessitates a partic- 
ular density of vortices, it can only occur when vortices 
proliferate, i.e. at the KT transition. This point is fur- 
ther verified by marking the percolation temperature on 
the curves for the PAD exponent, Fig. [5] For all values 
of disorder, this temperature corresponds to a ~ 0.38, 
midway between the values 0.25 and 0.5, exactly where 
the KT transition should occur. 

If these two transitions describe the same physics, then 
a curious question arises. For the KT transition the crit- 
ical divergence of the SC correlation length £ is expected 
to be £(T) - e s/VTc/(T c -T) where Tc the critical tem _ 

perature and B is some non-universal number. On the 
other hand, near the percolation transition one expects 
~ (p c — p)~ u where p is the link density (i.e. the 
concentration of links for which Da is above the thresh- 
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FIG. 4: The joint distribution of the local SC correlations and 
the local vortex density (normalized to unity at each vortex 
density). The strong correlation between the two values per- 
sists with increasing disorder, with the same functional form 
(solid line). 



old), p c its critical value, and v the correlation length 
critical exponent. Identifying the two transition, i.e. the 
critical temperature T c with the percolation temperature 
T p , and the two critical behaviors, implies a very specific 
dependence of the link density on temperature 



p(T) = Pc + Ae 



(4) 



The dependence of p on T is displayed in Fig. [SJ along 
with the solid lines, given by Eq.[?J with the resulting val- 
ues of B in the inset. The excellent agreement between 
the observed link density and the result of Eq. |4] gives 
further credence to the intimate connection between the 
KT transition and the percolation transition. This pecu- 
liar dependence of the link density on temperature indi- 
cates that as the system approaches the KT transition, 
the local values of Jij decrease rapidly, as the number of 
free vortices increases. Alternatively, when temperature 
is reduced through the critical temperature there is an 
avalanche of correlations through the system. Such be- 
havior cannot be captured, for example, by a disordered 
X-Y model, where the values of the local Josephson cou- 
plings are fixed and do not change. 

It would be interesting to investigate the critical be- 
havior as one increases the thickness of the sample. Here 
one may expect a crossover from the KT behavior to the 
mean-field BCS description. Whether percolation still 
plays a role and how it is related to the BCS transition 
is a question that we plan to explore in the future. 

This research has been funded by the ISF. 
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FIG. 5: Link density p as a function of temperature T. Large 
dots indicate T p and agree well with p = 0.5 as expected from 
percolation theory. Solid lines are fits according to equation 
[4] Inset: exponents (B deduced from the fit). 
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